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ABSTRACT 


Plasma temperatures and total radiant intensities were measured in air 
plasmas created in the reflected shock region of a 12- inch diameter^ arc- 
driven shock tube. The temperatures were required to explore the differ- 
ence between the total intensity measurements and theoretical predictions 
which were evaluated at the gasdynamic property values. The temperature 
measurements utilized the following spectroscopic techniques: (l) measuring 

the integrated intensity of certain atomic lines and multiplets^ (2) fitting 
profiles to measured line shapes, and ( 3 ) measuring the spectral intensity 
of the continuum at ^935 A. Temperatures between 10,750-13^000 K were measured 
with a typical precision of + h^o. This precision could not be maintained at 
higher temperatures. The results indicate that, on the average, the measured 
temperatures were below the gasdynamic values which could account for the 
difference between the total intensity measurements and the predictions in 
this temperature range. Further, the erratic shot to shot deviations of 0 
to 8^ between the measured and gasdynamic temperatures suggests that the actual 
plasma temperatures were not adequately predictable from the incident shock 
velocity. 
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FUJAL REPORT 
Contract MS 7-579 

RADIANT ENERGY TRANSFER MEASUREMENTS IN AIR 


INTRODUCTION 

The work done under this contract can be split into two distinct areas: one 

involves measurements and theoretical predictions of the spectrally- integrated 
radiant intensity of high- temperature air while the other concerns spectral 
measurements obtained using various spectroscopic techniques. Because the 
latest results of the spectrally- integrated area have already been published,^ 
this effort is briefly summarized below. The results and conclusions obtained 
serve as the introduction to the spectroscopic effort which occupies the bulk 
of this report. 

The measurements in both areas were performed in the reflected shock region 
of a 12- in. diam. shock tube equipped with an arc-heated driver. The test 
gas was air at = 0.2 torr; the incident shock speeds were 6. 5-9* 5 mm/iOisec 

and yielded reflected shock plasmas of 10 - 16^000 K, 2-5 atm, and p/Pq ^ 0.033 
Amagats. 

The spectrally- integrated effort involved measurements of the radiant intensity 
from these plasmas at paths of 3-12 cm using energy detectors in both a window- 
less and a window- type configuration. The salient results are shown in the 
two figures located in the appendix of this report. These show a comparison 
between the experimental results and theoretical predictions for both the 
total (i.e. windowless) radiant intensity (Fig. A-l) and for that fraction 
transmitted by a quartz window (Fig. A-2). Also included in the appendix are 
tables containing numerical values for the radiant intensity predictions 
(Table A-l) and the reflected shock Hugoniot states (Table A-2) on which the 
figures were based. A more complete description of the experimental apparatus 
and technique and the origin of the theoretical predictions is given in Ref, 1. 
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The comparisons between experiment and theory as presented by Figs. A-1 and 
A-2 show essentially the same behavior: good agreement at the lowest incident 

shock velocities with the experimental data falling steadily below the pre- 
dictions and reaching a factor of two below at the highest velocity. It is 
important to appreciate that the experiment measured the radiant intensity as 
a function of incident shock velocity while the predictions were evaluated at 
the thermodynamic states obtained from gasdynamic considerations. Therefore 
these comparisons between theory and experiment must include a consideration 
of the validity of the gasdynamic temperature. 

This basic question^ which precluded a meaningful comparison between theory 
and experiment, motivated the spectroscopic effort which, in fact, occupied 
the bulk of the past contractual period. This work had two important objectives. 
The first was to measure the plasma temperature directly and hence resolve the 
fundamental difficulty discussed above. The second was to generate spectral 
information about high- temperature air radiation to permit a more discerning 
evaluation of current radiation theories. The effort was concentrated on the 
former and more important objective. For reasons to be discussed later, the 
plasma temperatures were only measured to 13^500 K. However even these measure- 
ments revealed significant deviations between the actual and the gasdynamic 
plasma temperature and further work in the area is continuing. 

The temperature measurement problem can be placed in perspective by noting 
that the total intensity for a 5 cm path can be represented by la T^ 

over the conditions of interest. The temperature exponent is about 12 at 
10,000, 9 at 13^000 and 5*5 at l6,000 K. Thus a 30^ change in radiant inten- 
sity corresponds to a 2^ change in temperature at 10,000 K and a 5*^ change at 
l6,000 K. Thus it is clear that precise temperature measurements were re- 
quired. 

This high sensitivity of radiation to temperature, which requires precise 
temperature measurements, can be used to advantage by determining a temperature 
through measurements of radiative quantities. Two that proved suitable were 
the absolute intensities of certain well-known atomic lines and the shapes 
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(i.e. the half -widths) of these and other lines. As discussed later, a line 
intensity measurement yields a temperature based on a Boltzmann distribution 
of the upper state population while a shape measurement yields the electron 
density which is related to temperature through the Saha equation. 


The measurement of these two quantities required quite different instrumenta- 
tion and techniques. The line shapes were determined by photometric techniques 
via time- integrated spectrograms taken using film and a rapid-closing shutter 
in a spectrograph. Photographic recording captures a wealth of information 
about relative spectral intensities which is essential for determining line 
shapes, however it is not amenable to measuring absolute spectral intensities. 
These were measured over narrow spectral intervals using photoelectric record- 
ing techniques in a polychromator (i.e. a multi-channel monochromator) and 
were, of course, time- resolved measurements. The features of these two spectro- 
scopic techniques are simmiarized below: 


Instrument 


Recording 

Technique 


Time 

Variation 


Type of 
Measurement 


Spectral 

Variation 


Spectrograph 

Polychromator 


Photographic 

Photoelectric 


Time- inte grated 
Time- resolved 


Relative 

Absolute 


Spectral 

Spectrally- 

integrated 


As noted, the techniques complement each other and the combined results can 
yield a fairly complete picture of the plasma radiation. In fact, as will be 
explained later, the results of both techniques were required to obtain 
temperatures from both line intensity and line shape measurements. 


DESCRIPTION OF APPARATUS AND TECHNIQUE 

Both types of spectroscopic measurements were made on the same kind of instru- 
ment: Model 78-000 Jarrell-Ash 1.5 meter Wadsworth spectrographs. Each con- 

tained a 590 groove/ram grating and could cover from 2100 to 78 OO A at one 
time in first-order with a reciprocal linear dispersion of about 11.2 A/mm. 

The ruled grating areas were 5^ x 40 mm which means that the instruments were 
effectively f/29 and had a theoretical resolution of O.O 7 A at 2200 A and Q.2k 
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A at 7800 A. These parameters would dictate a spectroscopically-optimum 
entrance slit width of about 17 for recording at 6,000 A. In practice 
resolution was sacrificed to gain light and 50 entrance slits were used. 

Since these 50 entrance slits were from 2 to 8 times wider than the width 
for optimal resolution, diffraction effects were small and the resolution was 
very near the spectral width of the entrance slit (0.55 A) as will be dis- 
cussed shortly. One of these instruments was used as a conventional spectro- 
graph while the other was modified for use as a polychromator. The spectro- 
scopic techniques involved in these two different usages will be discussed 
separately. 

The terminology is however the same and begins with the spectral intensity 
(i^e. W/cm -sr-|ji) emitted by the plasma. Some of this radiation 

is collected by the spectrograph, dispersed, and brought to a spectral focus 
along the focal plane. Here it is characterized by the focal plane heat flux 
(i.e. q(^) in W/cm ) which varies with distance along the focal plane or, 

more conveniently, as wavelength. The role of an exit slit will be discussed 
later. A film placed here receives this heat flux over an exposure time with 
the result that a certain energy per unit area is deposited on the film. By 
the photographic development process, this energy per unit area is converted 
into an optical density which, by photometric techniques, is converted back 
into a spectral exposure (i.e. e(X) in erg/cm ) or a relative spectral 
exposure [e(x)/e(?\.^)] . At this point a distinction is made between relative 
spectral exposure and intensity . These are related through the instrument 
function which will be discussed later. First the methods by which the 
relative exposure was obtained from an exposed film will be discussed. 

Homochromatic Photometry 

The spectrograph was used to obtain relative spectral exposures by the technique 

2 

of homochromatic photometry (ref. Chap. 10 of Sawyer or Chap. 7 of Harrison 

o 

et al.*^). This technique essentially assumes that the wavelength intervals 
under consideration are sufficiently narrow that film sensitivity variations 
with wavelength are negligible. Otherwise, conventional photometric techniques 
are employed. 
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A photographic emulsion;^ upon being exposed and developed^ contains particles 
of silver which make the negative appear dark or even black when viewed by 
transmitted light. This darkness is characterized by an optical density 
which is the logarithm of the reciprocal transmittance. The relation between 
optical density and the amount of exposure is expressed as a density- log 
exposure curve commonly known as an H-D curve. Since optical densities are 
measured after the fact^ an H-D curve can only be obtained if the negative 
contains areas of known exposure and^ for relative work, these exposures 
need only be known relative to one another. The whole technique of homo- 
chromatic photometry involves obtaining and using the appropriate H-D curves. 
However, this is not trivial because the H-D curves depend on a great many 
parameters including the wavelength and the rate, amount, and time of exposure, 
the age, history and of course type of film, and the many variables (time, 
temperature, type, -concentration, etc.) of the development process. Good 
technique requires that the H-D curves be obtained under conditions as close 
as possible to those of the actual experiment. 


Because the Wadsworth mount is stigmatic, this was readily accomplished by 
placing a calibrated neutral density step wedge directly in front of the 
entrance slit during the exposure to the plasma radiation. The resultant bands 
of reduced exposure are shown by the following sketch of a typical negative. 



Neon reference lines 
Mercury reference line 
Unattenuated exposure 

Six attenuated steps 
Unattenuated exposure 


Sketch of a t3npical spectrogram showing the various 
exposure bands and the reference wavelength lines. 
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The step -wedge -was mounted on a special plate -which replaced the Hartmann 
diaphragm typically found on most spectrographs. The step wedge had six, 
metal- film attenuator steps on a quartz substrate. The steps had optical 
density differences of about 0.2 and occupied 8 mm of a 15 inm tall entrance 
slit. These dimensions are reflected in the above sketch because the spectro- 
graph had unity magnification. 

This sketch also indicates how the entire 15 nna image distance -was utilized 
to record an unattenuated exposure and reference wavelength information. The 
technique was to remove the dark slide from the spectrograph camera shortly 
before each shot. Then a mercury spectrum was put over the entire height of 
the entrance slit followed by a neon spectrum over the extreme edges. These 
spectra were obtained from low-pressure discharge tubes and the exposure 
times were several seconds. The purpose of the neon spectrum was to provide 
additional lines in the interesting red region where very few first-order 
mercury lines exist. Then the step wedge was inserted, the shock tube fired, 
and finally the dark slide -was replaced and the film removed and developed. 

By this technique, each spectrogram contained all the wavelength information 
needed to obtain the dispersion and absolute wavelengths and all the exposure 
information to obtain the H-D curve of the film at any desired wavelength and 
hence the relative exposure profile. This latter process involved several 
steps which are outlined below. 

As shown above, the six steps of the step wedge produced bands of reduced 

exposure across the film. The amount of reduction in each band depended on 

the optical density of each step and determining quantitative values can lead 

to much confusion. Optical density is defined by the logarithm of 

where I is the amount of light incident and I the amount transmitted, 
o 

The confusion lies in how the transmitted light is collected. Since some 
filters mainly scatter light (photographic emulsion), some absorb (gelatin), 
while others mainly reflect (metal film), the solid angle used to collect 
the transmitted light is important. Thus optical densities are further 
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qualified by the adjectives specular or diffuse. The correct density for 
the present application was that obtained by collecting the light transmitted 
into an f /29 solid angle - the same as the spectrograph. 

To accomplish this, the spectrograph entrance slit (with the step wedge in 
place) was illuminated by a tungsten strip lamp. A small slit was placed 
on the envelope of a photomultiplier tube and the slit- tube assembly was 
carefully located in the focal plane at various wavelengths. Then this 
assembly was racked up and down (ref. above sketch) across the various bands. 
The ratio of the outputs gave the transmittance of each step: they were 

0 . 63 ^ 0.40, 0 . 27 ^ 0 . 115 , and 0.077 and were uniform from 4500-7000 A. 

The optical densities of each spectrogram were measured using a Jarrell- Ash 
Model JA 2310 Microdensitometer equipped with a strip chart recorder. The 
microdensitometer output and hence the pen deflection was linearly pro- 
portional to the transmittance T of the emulsion which is related to 
density by D = 2 - log^Q(T) with T in percent. By using special density 
paper (Bristol Chart No. R0333) in the recorder, the algebraic conversion 
step was eliminated and density was read directly (albeit nonlinearly) from 
the paper. The microdensitometer slit was typically set at 0.7 3nm tall to 
avoid the edges of the exposure steps. Since the Wadsworth mount is only 
truly stigmatic at the grating normal (5IOO A), the edges became obscure at 
wavelengths far from the grating normal. The width was a compromise value of 
25 p- or roughly half the physical slit width of the spectrograph. A wider 
microdensitometer slit would integrate across more film grains and hence 
reduce noise, but this would also widen the instrument function which will 
be discussed shortly. The scanning rate was kept low (equivalent of 5*5 A/ 
min) to allow time for the recorder pen to follow all the spectral details. 

H-D curves were constructed by measuring the optical densities of the various 
exposure bands at the wavelengths of spectral features of interest. The basic 
premise of homochromatic photometry is tantamount to assuming that each H-D 
curve applies over, say, the full width of the line or multiplet under con- 
sideration. Therefore, the microdensitometer scan along any exposure band 
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could be converted^ through the H-D curve, to an ordinate of relative exposure. 
Normally the unattenuated exposure band was used for this procedure since the 
negatives were generally underexposed and it is desirable to work on the most 
linear portion of the H-D curve. The reference wavelengths included on the 
microdensitometer scans allowed the wavelength at any point to be determined 
by linear interpolation. In practice the wavelengths were often read directly 
from a preset Gerber scale. Thus the final result of this technique was a 
plot of relative exposure versus wavelength - a relative exposure profile. 
Examples are given later in this report. 

Instrument Function 

The relation between the exposure profile discussed above and the spectral 
intensity profile emitted by the plasma is described by the instrument function. 
The spectral features of the emergent intensity were broadened both in the 
spectrograph itself since at best the focal plane contains a series of over- 
lapping monochromatic images of the entrance slit and in the microdensito- 
metric scanning process where again another slit of finite width was employed. 

Suppose that the entrance slit (area A ) was uniformly irradiated by a steady 

intensity l(^ ) from a monochromatic source at wavelength \ . Assuming 

that the grating represents the limiting aperture, the spectrograph subtends 

a solid angle A/f where A is the grating area and f is the slit-to- 

grating distance. ' Thus the product l(X ) A A/f represents the rate at 

s 

which monochromatic energy enters the spectrograph and of this a fraction 

■X* 

t(X ) reaches the focal plane in first-order. This energy is distributed 
along the focal plane in a yet unknown manner, however the integral of the 
focal plane heat flux q(X) must equal the incoming rate, hence 

i(x*) t(x*) ^ f 

■'’The presence of collimating lenses or mirrors complicates this definition but 
should not cause difficulty by noting that the expression represents the solid 
angle that the instrument subtends at the entrance slit. 

8 



I 


where h is the height of the focal plane ^ ELD the reciprocal linear 
dispersion, and AX the wavelength increment over which q(X) is nonzero. 

(it is convenient to express lengths along the focal plane in terms of wave- 
lengths by the expression dX = ELD dx), 

Now if a spectrogram were taken of this monochromatic irradiation, the result- 
ant relative exposure profile is defined to be the instrument function, 

y •X’ 

p(X-X ) = e(X)/e(X ), where e(X) is the spectral exposure at wavelength 
X and is normalized to the maximum value e(X*). The instrument function 
for this experiment is shown by Fig. 1. It was measured at three wavelengths 
using the reference mercury and neon spectra along the edge of an actual shot 
spectrogram together with H-D curves obtained from nearby spectral features. 

The particular lines were selected because they were weak and the resultant 
film densities were comparable to that of the plasma spectrum- Bright refer- 
ence lines could easily become grossly overexposed and lead to broadening by 
film halation. It is noted that the values at each of the three wavelengths 
are quite similar and hence the instrument function was not dependent on wave- 
length. Close scrutiny shows that the broadest function was obtained from 
the reddest line which should be expected since diffraction effects are least 
negligible at long wavelengths. 

Since diffraction effects should be small because of the relatively wide 
entrance slit employed, the instrument function would be expected to resemble 
the result of scanning a geometric image of the 50 li- entrance slit with a 
25 Pi microdensitometer slit; i.e. a triangular function with a base of 100 |Jb 
or 1.12 A. Since the actual function was trapezoidal with a slightly narrower 
base, the actual microdensitometer slit was probably slightly smaller 22 p^). 
More important is the elimination of other types of broadening phenomena such 
as improper focus, optical aberrations, etc. This predictability of the 
instrument function means that the monochromatic energy was uniformly distri- 
buted over the geometric image of the entrance slit in the focal plane. This 
will be important in the polychromator discussion which will follow shortly. 
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The equation relating the absolute exposure at a wavelength X resulting 
from the steady monochromatic irradiation of the entrance slit at a wavelength 
A.* is obtained from Eq. (l) by introducing an exposure time t and the 
definition of the instrument function with the result 

e(x) , T(X*)t 

*' y 

A\ 

For a stigma-tic instrument with unity magnification, the factor A /h becomes 

s 

the physical width of the entrance slit. 


The extension of this concept to non- mono chromatic intensities is straight- 
forward. Any spectral intensity broken into a series of 

"monochromatic" lines by multiplying by a differential wavelength increment 
(±\ and treating each incremental line in the above manner. The resultant 
exposure at any wavelength is the sum of all possible contributions and is 
given by 

a(X) T(X)t ^ 


*AX 


( 3 ) 


Here the transmittance t(X) was assumed to remain constant over the instru- 
ment function. This equation is generalized by invoking the well-known 
principle that the intensity of a source cannot be increased by optics. 

Thus the intensity in the above equation is that emitted by a source provided 
that (l) the transmittance includes that of the optical system and ( 2 ) the 
grating and the slit remain the limiting apertures. 


The above equation relates the spectral intensity emergent from the plasma 
to the absolute spectral exposure. The techniques of homochromatic photo- 
metry yield relative spectral exposure profiles. These are expressed in 
terms of the instrument function by normalizing the above equation with 
respect to a peak exposure e(X^). Since it was shown that the instrument 
function was independent of the absolute wavelength, Eq. (3) becomes: 
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e(X) ^ AX 

f I^(X») F(X^-X»)dX» (^) 

*AX 

If changes in the spectral intensity over the vidth of the instrument function 
are neglected^ this equation reduces to 

e(X) -V’') 

which is to say that the relative exposure profile is equal to the relative 
spectral intensity profile provided that the instrument function is negligibly 
narrow. 

In the preceding discussion it was tacitly assumed that the plasma intensity 
remained constant throughout the exposure time. The constancy was verified 
on each shot by a photoelectric monitor of the zero- order light in the spectro- 
graph and also by the simultaneous signals from the various polychromator 
channels. Examples of these traces will be presented later. 

Spectrograph Optical System and Rapid-Closing Shu tter 

The spectrograph was coupled to the shock tube by the optical system shown 
in Fig. 2. This optical system was designed with two requirements in mind: 
firsts the accepted beam had to be reasonably narrow over the full 30 dia- 
meter of the shock tube to avoid viewing the end wall and interface. Secondly, 
the entrance slit (or, more precisely, the step wedge) had to be uniformly 
irradiated so that the technique of homochromatic photometry would be valid. 
This latter requirement was met by using the 10 cm f.l. lens (ref. Fig. 2) to 
focus an image of the entrance slit on the l4 cm f.l. lens. Since this latter 
lens was uniformly irradiated by the reflected shock plasma, then so was the 
entrance slit. The accepted beam was made narrow by employing the cylindrical 
lens at slightly more than a focal length away from the 1^ cm f.l. lens. The 
centerline of the beam was located 7 from the end wall. 
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Fig. 2 Optical System for Spectrograph 



This optical system used ordinary fused quartz lenses throughout and trans- 
mitted down to at least 2500 A (this was probably the film cut-off caused 
by gelatin absorption). It was not achromatic and hence was focused in the 
green 5500 A) as a compromise over the visible region- It did go out of 
focus in the ultraviolet but this was of little consequence since this region 
was used for survey purposes only. The field and aperture stops were placed 
at images of the grating and slit to prevent excess light from entering the 
spectrograph and causing possible scattered light problems. These stops 
were slightly larger than the size of the visible images so that they were 
not limiting apertures. The spectrograph was fully illuminated in the visible 
region by this optical system. 

In practice the spectrograph was left open from several seconds prior to each 
shot until about 30 l-isec after the reflected shock passed through the field of 
view. The exposure from the ambient light and the incident shock was negligible 
in comparison to that from the reflected shock plasma- Thus the shuttering re- 
quirements were for a rapid-closing^ capping-type shutter. 

U 

The concept of this shutter was obtained from the article by Wurster although 
the final design was somewhat different. Basically the shock wave and/or debris 
from an exploding wire are used to bend a flap of metal foil through 90 ^ and 
hence cover an aperture. Figure 3 contains several sketches to show how the 
wire was mounted and how the effects of the explosion were directed toward the 
flap. The entire assembly shown was mounted in a brass case 75 diameter by 
to mm thick. The foil was ordinary "Super Strength Alcoa Wrap" which had a 
thickness of 0-7 mil. The 3 niil diameter by 15 mm long copper wire was 
exploded by switching a 4 iJifd capacitor charged to about 3*3 kV via a 5C22 
thyratron tube. The shutter would typically close in 5 M-sec with closure 
occurring about 35 M'Sec after a trigger pulse was sent to the thyratron grid. 

The initial firing of the thyratron caused an electromagnetic disturbance in 
the other instrumentation. Therefore the timing was such that the thyratron 
was triggered before the reflected shock reached the field of view of this 
instrumentation and hence this disturbance occurred and was over before any 
reflected shock information was recorded. 
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The action of the shutter was recorded for each shot hy the photoelectric 
monitor of the spectrograph zero- order light. Thus the exposure time and 
the temporal behavior of the radiation during the exposure were known. 

Examples of this monitor will be given later. No radiation emitted by the 
shutter was discernable on any spectrogram and this is probably explained by 
the physical location of the shutter components; the exploding wire is sur- 
rounded by the wire cavity until the flap opens and hence the shutter closes. 

All the spectroscopic work was done with a first-order spectrum. The higher- 
order spectra were blocked by placing a strip of Rohm and Hass Type II UVA 
Plexiglas along the focal plane from 38 OO A to beyond 7 OOO A. This material 
is opaque below 3365 A and transmits well above 3^60 A. Thus the combination 
of quartz lenses and Plexiglas cleared the first-order spectrum to about 68 OO 
A. A Wratten #8 (K2) filter was used occasionally to clear to even longer 
wavelengths^ but this appeared to cause interference problems when used near 
the focal plane and^ since the film sensitivity dropped sharply above about 
6900 Ay this filter was seldom used. 

Spectrograph Film and Exposure 

The film used throughout the experiment was Kodak 2^75 Recording Film. This 
is an extremely high speed (ASA 1000 ) panchromatic film with extended red 
sensitivity. It was developed for 5 minutes in room- temperature DK -50 follow- 
ing the manufacturer’s recommendation. Unfortunately high speed film is grainy 
and 2^75 film is rated as a coarse-grained film. Attempts to use a slower, 
finer-grained film (Kodak 2^96 RAR film) met with no success because of in- 
sufficient light. 

Although the spectrograph was used to obtain relative rather than absolute 
exposures, inevitably one wishes to work on an approximately absolute basis 
if only to consider exposure details. The necessary equation is obtained 
from Eq. ( 3 ) by neglecting the effect of the instrument function with the 
result 
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e(X) 


w(rld)a 


t(\) 


(5) 


For the JACO !• 5 spectrograph with a. ^0 slit, the hracketed factor is 
- R o 

5.U X lO" jJL-sr (a = UO X 54 imn , W = 50 p-, RID = 11-2 A/mm, f = 1.5 m)- 


The transmission of the spectrograph was obtained by illuminating the entrance 
slit with monochromatic light obtained from a monochromator through a system 
of lenses and apertures with a larger f /number (i.e- a smaller ‘solid angle) 
than that of the spectrograph. Thus the ratio of the light reaching the 
focal plane to that measured behind the entrance slit gave the transmission 
(for first-order) of the spectrograph. This was measured with a linear photo- 
multiplier tube and the result was that the transmission reached a maximum 
value of 60^o at 3500 A and dropped steadily to about 30^ at 7000 A. The 
transmission of the optical system was determined by back- lighting the spectro- 
graph and determining the ratio of the light emergent from the last lens to 
that emerging from the entrance slit. This ratio was about 50^ at all wave- 
lengths. The product of these two numbers thus yields the transmission re- 
q^uired in the above equation. 


To measure the sensitivity of the Kodak 2475 Recording Film, a 0.l6 sec 
exposure was taken of a calibrated tungsten strip lamp - a source of known 
intensity. The resultant spectrogram was densitometered at a number of wave- 
lengths to obtain the H-D curves which, in this case, were used to yield the 
relative exposure necessary to obtain a density of 0.3 above fog. The absolute 

value of the 100^ exposure was calculated from the above equation using the 

known lamp intensity and the measured transmission. The results are shown on 
Fig. 4 together with the Kodak curves for this and two other films. As can 

be seen, the data points agree quite well with the manufacturer's information 

and a log sensitivity of 2 is typical for 2475 film. Whether this value of 
film sensitivity is valid for a calculation of an intensity from the re- 
flected shock plasma at a typical exposure time of 30 pisec depends on the 
reciprocity behavior of the film. The reciprocity law states that the film 
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Fig. h Spectral Sensitivity of Various Kodak Filins 


sensitivity is dependent only on the product of the exposure rate and time 
and is independent of either factor taken separately. Milne and Eyer^ could 
find no reciprocity failure for Kodak Royal-X Pan Recording Film (a type 

-Ij. -8 

similar to 2475) Q*t exposures from 10 to 10 sec so the neglect of 
reciprocity failure for this calculation appears to he valid. 

—2 2 

Thus ve assume that an exposure of lO" erg/cm is required for a density of 
0.3 above fog (a minimal exposure). Using the above equation and assuming a 
transmittance of 20^^ the intensity required for a 30 ijisec exposure is I = 3100 
W/cm -sr-jJi. For visible radiation and a 30 cm pathlength in the reflected 
shock plasmas of this experiment, this spectral intensity is typical of the 
radiation at line centers from an 11,000 K plasma and continuum radiation 
from plasmas above 12,500 K. 

Thus it is clear that the experiment did not enjoy excess light and the finer- 
grained but less- sensitive (by factor of lO) Kodak 2496 Film could not be 
employed. Recently another extremely fast Kodak film has become available - 
Kodak 2485 High Speed Recording Film (formerly Type SO-I 66 ) - and promises to 
be nearly 10 times faster than 2475 at 4000 A. It -was not used in this experi- 
ment, although it will probably be tried in the future. In an attempt to 
capture the near- infrared region, a short length of Kodak High Speed Infrared 
Film was placed in the 7 OOO- 78 OO A portion of the spectrograph camera. This 
met with no success (negligible density change) and the cause appears to be 
reciprocity failure in the red- sensitive emulsion since good exposures were 
obtained from the standard lamp. However quantitative work needs to be done 
before definite statements can be made. 

Polychromator 

The polychromator was used to measure photoelectrically the absolute intensi- 
ties averaged over 5-50 A spectral intervals at selected wavelengths in the 
visible spectrum. As mentioned earlier, the basic instrument was another 
JACO 1.5 m Wadsworth Spectrograph. A large (l x I .5 x 2 ft), light-tight box 
was added at the camera end so that phototubes could be conveniently located 

Kodak literature for the 2475 film (Pamphlet No. P-95) indicates that 
deviations between 10"'^ sec and O .16 sec are negligibly small (^ 10 ^). 
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behind the focal plane. The polychromator was coupled to the shock tube by 
the optical system shown on Fig. 5* Again the requirement was for a narrow 
beam along the tube diameter but in addition the system had to form a small 
image of the entrance slit in the center of the tube for calibration purposes. 
The overall length of the system was much longer than that of the spectrograph 
(Fig. 2 ) because of physical space limitations. Glass lenses were used because 
radiation below 3500 A was of no interest here. The polychromator field of 
view was directly opposite to that of the spectrograph so that both instru- 
ments viewed the same gas sample^ but from different ends. 


F±-vey fixed- type exit slits were mounted along the focal plane. So that the 
photocathode of each photomultiplier tube would be uniformly irradiated regard- 
less of the spectral energy distribution across each exit slit^ the system 
sketched below was employed. 



23 mm diam 
glass achromat 


Top view of component arrangement behind focal 
plane of polychromator , 

Here light emergent from each slit was collected by the 17 cm f.l. achromats 
which foimied l/8-size images of the grating on the photocathodes of 1 P 28 
photomultiplier tubes. Corning glass color filters were used to block 
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TOP VIEW 


Fig. 5 Optical System for Polychromator 



scattered light and (in one instance) second- order spectra. The open spaces 
between each exit slit were covered with black photographic tape so that the 
only light entering the light-tight box had to come through an exit slit. 

The phototubes were collimated as shown to eliminate possible cross talk. 

The phototubes were culled from a large supply by selecting those which 
exhibited the best signal/noise ratios at the wavelengths of interest. 

The factors governing the selection of the center wavelength and spectral 
bandpass for each exit slit will be discussed later and here these values 
are assumed to be given. The technique used to locate these slits was to 
have a second alignment slit cut in the same piece of metal. A nearby^ 
strong reference line from the spectra of Hg^ ky He^ or Ne was selected for 
this alignment. This line^ together with the above values and the measured 
(photographically) reciprocal linear dispersion gave the necessary input to 
design an exit slit similar to that sketched below. 



Sketch of fixed exit slit used in polychromator 


The exit slits were milled from 0.010- inch hard aluminum sheets about I .5 inches 
square. The exact widths and separation were determined by measurement on an 
optical comparator. 

The slits were clamped along the focal plane by screws and washers to permit 
adjustment. Alignment was accomplished by locating the appropriate reference 
line in the field of a short- focus telescope and then manipulating the exit 
slit until this reference line was in the center of the alignment slit when 
the exit slit was locked tight. With a 50 M* entrance slit^ the width of this 
reference line was 50 iJi or 2 mil. An easily observable misalignment would be 
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2 mil which is equivalent to 0. 5 A in wavelength along the focal plane. Thus 
the exit slits could be located in the above manner to better than + 0.5 A. 

Polychromator Calibration and Analysis 

The polychromator system was calibrated by replacing the unknown reflected 
shock plasma with a source of known spectral intensity. This was an Eppley 
Standard of Spectral Radiance - a General Electric 30A/T24/7 tungsten ribbon 
filament lamp calibrated by Eppley from 0.25 "to O .75 at 35 amps. In 
practice a similar lamp was used as a working standard. The lamps were 
placed at the center of the shock tube where the slit image was comfortably 
smaller than the filament. Thus the system was calibrated through the same 
lenses^ windows^ etc.^ as used for each shot. This did mean that the cali- 
bration had to be performed hours before each shot - before the tube was 
closed and evacuated. However with one exception, the calibrations did 
not change significantly from day to day so this was not a serious disadvantage. 

To provide a periodic calibration light signal, a 36 OO rpm chopper was placed 
immediately in front of the entrance slit. The chopper blade was 5*2 inch 
diameter and contained a single slit O.O 5 O inch wide. This yielded a flat- 
topped 50 \isec light pulse every I /60 sec and hence provided a fair simula- 
tion of the time behavior of the transient light pulse from the reflected 
shock region. 

The amplitudes of the plasma light pulses were from ^0-^000 times greater 
than that of the standard lamp and thus the linearity of the photomultiplier 
tubes was an important consideration. Linearity was established over the 
full range by replacing the standard lamp with a Xenon flash tube and insert- 
ing various neutral density filters into the optical path. The amplitude 
disparity was slightly reduced by opening the polychromator entrance slit 
to 300 iJi during calibration (the slit image remained sufficiently narrow) 
which gave a sixfold increase in the energy received by each phototube. The 
principal advantage of this was to improve the signal/noise ratio of the cali- 
bration signals. The exact factor was measured by using the zero-order light 
and it appeared to hold for all the phototubes. 
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III! 


The equation for the rate at which energy passes through an exit slit on the 
polychromator focal plane is readily derived hy rewriting Eq. ( 3 ) in terms of 
a focal plane heat flux. However for this application the Instrument function 
describes the appearance of a monochromatically- Irradiated entrance slit in 
the focal plane. From the earlier results, this is merely the geometric image 
of the entrance slit and thus the focal plane heat flux becomes 


q(^) 


= 4 T(X) 

f 



( 6 ) 


This is the heat flux at a particular point on the focal plane and it receives 
contributions from wavelengths \ ± W*/2 where W’ is the spectral width of 
the entrance slit. 


The contribution from each wavelength is obtained by differentiating the above 
equation whereupon 


d-q.(^) = ^ t(x) i^(x)dx 

Here dq(X) represents the incremental heat flux resulting from the mono- 
chromatic intensity I^(X)dA.. This flux is unifonn over the entrance slit 
image which occupies a distance W or spectral interval W* = ¥(RLD) along 
the focal plane. The distribution of these differential heat fluxes is best 
shown graphically. The sketch on the next page indicates how an edge of the 
exit slit masks a portion of the images resulting from wavelengths near that 
of the exit slit edge. The situation for the other edge would^ of course^ 
be similar. This sketch considers intervals d\ = clarity. 
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If the edge of an 
exit slit were 
placed here: 
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Only these fluxes 
would pass 


Schematic representation showing passage of 
differential fluxes by an exit slit edge. 


To obtain the rate at which energy passes through an exit slit, the areas 
have to be modified to account for this masking effect. This is accomplished 
by the weighting function R(X) sketched below. 



Exit slit weighting function. 

Now the differential rate for the energy passed by an exit slit is given by 

de(X) = dq(X) A R(X) and, by integrating over the nonzero limits of R(X), 
s 

the total rate E is 
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(8) 


AA 

E(AX) = 



R(X) x(X)dX 


ax+w 


The weighting function only involves the intensities near the exit slit edges 
and these edges were placed in spectral regions where the variations of the 
spectral intensity with wavelength were small. By neglecting these small varia- 
tions^ the weighting function B(X) can be "squared up" and further by assum- 
ing the transmission to be constant over the entire exit slit widths the above 
equation can be simplified to yield 


vf 


e(ax) 


AA r 


I^(X)dX 


(9) 


x-fi 

o 2 


During the polychromator calibration the source was a standard lamp whose 
spectral intensity varied slowly and smoothly with wavelength. By choosing 
a mean value the following operational equations for the polychromator can 
be derived from Eq. ( 9 ): 


I,, = I, (std. = s y - at y (xo) 

AX AX 

Here I^(std. lamp) is the mean spectral intensity of the standard lamp and 

^plasma^^std lamp ratio of the shot signal to the calibration signal 

obtained from a given photomultiplier tube. Thus is the mean spectral 

intensity of the plasma averaged over the spectral width of the exit slit as 
given by the third equality. The fourth equality was obtained by inserting a 
normalizing constant to make the integrand resemble the relative spectral 

intensity profile as obtained by the spectrographic techniques discussed 
earlier. The above equations represent the limit of the simplification of the 
polychromator results without further knowledge concerning the spectral nature 
of the plasma radiation. 
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Combination of Eesults from Spectrograph and Polychromator 


Tlie spectrograplaic results provide this detailed spectral information and 
can he used to perform the integration required hy the last equality in Eq. 
(lO), For the moment, changes in the spectral intensity over the- width of 
the instrument function will he neglected so that the relative exposure 
profile can he considered equivalent to the relative spectral intensity 
profile. Thus the position of the exit slits can he drawn on the profile 
and, hy planimetry, an average relative spectral intensity (r) over AX can 
he obtained as sketched helow. 



the average relative spectral intensity. 
This operation yields, from the spectro graphic result; 



and combining this with the appropriate polychromatic result from Eq. (lO) 
yields 

r 

This is an important equation because it provides absolute calibrations for 
the relative spectral intensity profiles from the spectrograph without having 
to perform quantitative photometry. 
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PIASMA THEEMOMETRY 


In this chapter the equations needed to calculate a plasma temperature from 
the polychromatic and spectro graphic results are derived together with the 
appropriate error analyses. 

General Error Analysis 

Basically the plasma temperature was calculated from the measured value of a 
temperature- dependent radiative quantity M. Such a quantity^ emitted hy a 
steady isothermal plasma, will generally depend on the two thermodynamic 
variables which define the state of the plasma together with certain physical 
constants. By selecting temperature and pressure as the thermodynamic varia- 
bles and denoting the lumped physical constants as C, this functional rela- 
tionship can be expressed by M = f(T,P,C). The total differential of this 
relation can be written in the following nondimensional form: 

m Vm\ ^ C SMi ^ 

« M ' P - 

T T (13) 

« P,C 

Here the left side represents the uncertainty in the calculated temperature 
while the numerator on the right contains terms which describe the uncertainty 
in the measurement of M itself together with the effect of the uncertainty 
in the pressure and the lumped physical constants. The denominator will be 
called the temperature sensitivity and clearly it will mitigate the other 
uncertainties if it is large. 

Thus the selection of the quantities M used to determine the plasma tempera- 
ture depended in part on the temperature sensitivity. This in turn is a 
function of temperature and can limit the region of usefulness as will be 
demonstrated. For the lower plasma temperatures the shape and/or the total 
intensity radiated by certain atomic lines and multiplets were selected as 
suitable quantities. 
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Derivatives of Plasma Properties - 

To evaluate the temperature uncertainty resulting from the measurement of 
these quantities, certain partial derivatives of the plasma properties will 
he required. 


The Saha equation and the equation of state for the singly- ionized region of 
a neutral, monatomic (or fully dissociated polyatomic) plasma of a single 
specie are: 


2 U_L 

O' + 

1 - O' ~ U 
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A exp(- 1^) 


P = (2n^ + n^) kT 
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(14) 

(15) 


where a. is the degree of ionization, U^/U is the ratio of the ion to 
neutral electronic partition 'functions, I is the ionization energy, 

A = 4.83 X 10^^ cm~^ and n^ and n^ are the electron and neutral 

particle densities. These equations are given in Chap. 3 of Zel’dovich and 

g 

Eaizer. 


The requisite partial derivatives are obtained from these equations and can 
be expressed as follows: 
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Since 
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O' = f(T,P), these derivatives 


= 1 + 0' (18) 

are functions of temperature and pressure. 
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While the Saha equation is not strictly valid for air, it is noted that if 
U_^/U = 0«77 3 -nd I = 13.8 eV are used in Eqs. (l4) and (l5)^ the resultant 
values for the degree of ionization are within 1 ^ of the values calculated 
from the data in Table A2 where a = 0.01 to 0.5* This is excellent agreement 
for the intended purpose. Thus the above partial derivatives were evaluated 
along the path of the reflected shock Hugoniot states (ref. Table A2) by 
using I = 13-8 eV and the tabular values of the degree of ionization. The 
resulting curves are shown in Fig. 6 and will be used in the discussions 
which follow. The curves apply to either nitrogen or oxygen atoms because 
of the assumptions involved. 


Radiant Energy Transport Equations 

The basic equation for radiation transport in an isothermal plasma is 






(19) 


where is the specific intensity, 6 the pathlength, and the spectral 

absorption coefficient which contains additive contributions from all line 
L C 

and continuum (m<^) radiation. is the Planck function which, for 

wavelengths in microns, is given by 


^ 1.191 X 10 
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where = 1*^39 ^ 10 


( 20 ) 


Invearbing the basic transport equation (Eq. 19 ) yields 

I, 


rn L ^ C 1 , 
^^X = ^^"X ^X = ■ 6 
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( 21 ) 


For line radiation, the integrated absorption is defined as 
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( 22 ) 
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VALUE OF PARTIAL DERIVATIVES (DIMENSIONLESS) 


8 



Fig. 6 Partial Derviatives of Plasma Properties Evaluated Along the 
Reflected Shock Hugoniot 



where now refers to contributions from a particular line or multiplet 

only. For wavelengths in microns, the Integrated absoipition is: 
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( 23 ) 


where is the number density of the emitting specie in part./cm*^^ 

the partition function^ g the lower state degeneracy^ f the f-number^ 

and e- the lower state energy. The last three parameters were taken from 
^ 7 

the KBS tables whenever possible. These equations^ together with the com- 
bined results from the polychromator and spectrograph^ were used to determine 
a temperature from the measured intensity of a particular line or multiplet 
as outlined below. 


Thermometry of Line Intensity Measurements 

The required numerical work was done by a digital computer code originally 
written to analyze line radiation that was not optically thin. However since 
it is valid for any opacity^ its use became routine despite the fact that 
simpler expressions can be written when opacity effects are negligible. 

Basically the integrated absorption (Eq. 22) was computed from experimental 
measurements of a given line or multiplet. Then this value was used to 
determine a temperature by numerical inversion of Eq. (23)* 

The pertinent experimental results were the relative intensity profile from 
the spectrograph (see sketch on p. 27 ) together with the value of the normal- 
izing intensity as obtained from the combined result from the spectrograph 
and polychromator (Eq. 12). This experimental profile contains contributions 
not only from the line or multiplet under consideration but also from the 
continuum and possibly other multiplets. These contributions must be separated 
and this becomes increasingly difficult at the higher temperatures. Examples 
of this separation will be given in the next chapter. 


The continuum^ as used here^ means an unchanging spectral intensity and 
could very well include the underlying wings of other multiplets. 
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For the moment, consider this separation accomplished and select two wave- 
lengths and X^ on opposite sides of a given multiplet where the 

spectral absorption coefficient of the multiplet is effectively zero and 

Q 

where the continuum Intensity (ij^) can he determined. Equations (2l) and 
(22) can now be combined to yield 




(^)(^) 


dX - (x^-x^) 


where q 

- j « (. - K^,] 

Here spectral changes in the Planck function, continuum intensity, 
normalizing factor have been neglected. 


and 


(24) 


The computational scheme was to input the normalizing factor (Eq. 12), the 
relative intensity profile, and a first estimate of the plasma temperature 
which was used to evaluate the Planck function. The computer calculated a 
value for the integrated absorption (Eq. 24) utilizing Simpson* s rule and 
then, by iteration, found the temperature that satisfied Eq, (23)- This 
temperature was generally not equal to the first estimate, so a second 
estimate was made and this procedure continued until the two temperatures 
agreed. The printout included the continuum spectral absorption coefficient 
(|Ji^) and the peak emissivity these are noteworthy because of the 

consistency checks they provide as discussed later. These computations were 
performed on a time- sharing digital computer (TYMSHAPE). 


Up to this point the effects of the instrument function and the exit slit 
weighting function have been tacitly neglected. The exit slit weighting 
function can be ” squared up" (Eq. 9) without introducing a significant error 
provided the spectral intensity did not change rapidly at the exit slit edges. 
The slits were located so that this was generally true. The effect of the 
instrument function can be assessed by first recalling that broadening does 
not destroy energy. Thus in general the main effect of treating a relative 
exposure profile as a relative intensity profile is that the normalizing 
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intensity (Eq. 12) is too low. This error is partially cancelled in 

the subsequent spectral integration process and it can he shown that there 
is no resultant error when the radiation is optically- thin. When opacity 
effects become important^ the errors do not completely cancel but at the same 
time, opacity effects are concurrent with higher temperatures, electron densi- 
ties, and hence broader lines whereupon the instrument function itself becomes 
less important. Thus in summary the effects of the instrument function and 
exit slit weighting function were negligible when spectrally- integrated 
quantities were used as was the case in this thermometric technique. 


Error Analysis for Line Intensity Thermometry 

To obtain analytic expressions for the temperature uncertainties, the radiation 
must be assumed to be either optically- thin or optically- thick. Since the 
experimental conditions were more nearly optically- thin and since it can be 
shown that the optically- thick case yields only slightly greater temperature 
errors, only the optically- thin case will be treated here. 


The optically-thin assumption (l^ « or |i ^6 « l) permits the integration 
of Eq. (24) by the expansion In(l-x) = - x with the result: 


r 

J 


(25) 


The measured quantity M is considered to be the integrated line intensity 
The continuum contribution will be treated as a constant because its 
magnitude was not used to determine a temperature in this thermometric 
technique. 


Combining Eq. ( 23 ) with the Planck function (Eq. 20) yields the integrated 
line intensity 


= 1.056 X 10 
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(26) 



where is the upper state energy. Substituting this expression into Eq. 

( 13 ) and performing the required differentiation yields the following expres- 
sion for the uncertainty in the measured temperature : 
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(27) 


The temperature sensitivity of this technique is given by the denominator 
which represents the change of the excited state population with temperature. 
For most visible transitions of NI or 01^ the upper state energy is 

about 13 eV. The derivative of the partition function is small (?« O. 5 ) and 
changes slowly. However the first term^ representing the change in the 
neutral number density, is negative and becomes large at the higher tempera- 
tures (ref. Fig. 6). Thus the temperature sensitivity decreases rapidly at 
the higher temperatures and the result of this will be demonstrated shortly. 


The first term in the numerator describes the uncertainty in the experimental 
measurement of the integrated line intensity which, for optically- thin con- 
ditions, is given by Eq. (25)- However the actual determination is better 
described by referring again to the sketch on p. 27- The polychromatic 
result provides the absolute value of the integral of the spectral intensity 
over the width of the exit slit. From Eq. (lO), the uncertainties in this 
value stem from those of the standard lamp intensity (± 2^), the measurement 
of the standard lamp signal voltage (+ 3^)^ and the measurement of the plasma 
signal voltage (+ 5 ^) and when added together, give a total uncertainty of 
± 10^0 to the polychromatic result. 


Then two corrections were applied to this result and both were derived from 
the relative spectral intensity profile. First an additive correction was 
applied to account for the line wings cut off by the exit slit and, for the 
wide slits used herein, this was a relatively minor correction. The second 
and more important correction was a subtraction of the continuum contribution 
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and this involved first an estimate of the apparent continuum level from the 
relative intensity profile. Through experience gained during the shape- 
fitting procedures discussed next, this level could he picked to + 0.02 units 
of relative intensity. Finally deviations between this relative intensity 
level and the actual plasma emission must be considered. By definition this 
deviation is zero at = 1.0 but lower values involve the H-D curve 

and the attendant effects of film noise, weak exposures, and the assumptions 
of homo chromatic photometry. Quantitatively, these effects are described by 
the following equation derived from practical experience : 


A (^) = ± 0.04 (1 - 


(28) 


\o 


"Xo 


where the left hand side represents the deviations between the actual plasma 
emission and the relative intensity profile. 

The effect of these two corrections was quantified by assuming a reasonable 
line shape and exit slit width and various continuum levels. The results 
indicated that these spectrographically- derived corrections to the polychromatic 
result added another + 10^ uncertainty if the continuum level was below 0-5 

and that the uncertainty rose rapidly if the continuum level was higher. 
This latter behavior occurs because the process of subtracting the continuum 
yields a result for the integrated line intensity which is the difference 
between two large numbers, each of which has an associated uncertainty. 

Thus the first term in the numerator of Eq. (27) is estimated to be + 20^ at 

low continuum levels. The second term contains a derivative which is near 

unity (ref. Fig. 6) and the pressure itself which is considered uncertain to 

+ 10% since it was not measured during this experiment. The third term 

7 

represents the uncertainty in the f-number which for KBS grade B lines is 
± 10%. Hence the numerator of Eq. (27) becomes ± h0% and, when divided by 
a typical temperature sensitivity of 10 at 12,000 K, indicates a temperature 
uncertainty of ± k%. 
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At 15;>000 the teiriperature sensitivity drops to about 4*5 vhich would yield 
a temperature uncertainty of + 9^ even if the continuum level stayed below 
50 *^ which^ as will be shown later, it does not. Hence the thermometric 
technique of measuring the integrated intensity of an HI or 01 line is 
fairly precise for temperatures up to about 13,000 K, but this precision 
is rapidly lost at higher temperatures both because of a loss in temperature 
sensitivity and because of the high continuum levels. 

Thermometry of Line- Shape Fitting 

The shapes of well-known atomic lines provided the other plasma thermometer 
used during this experiment. The fitting process was done essentially by 
choosing a set of plasma properties, computing the resultant relative intensity 
profile, broadening this by the instrument function, and comparing this against 
the measured relative exposure profile. Then the plasma properties were varied 
until a best fit to the experimental shape was obtained. This procedure in- 
cluded the effects of finite opacity and involved the whole experimental shape 
rather than just the half-width. The derivation below applies to a single 

■X* 

line from neutral atomic nitrogen which has a Lorentzian shape. Other atoms 
and shapes were dealt with in a similar manner. 

The line shape is normally characterized by a half half-width 7 . For neutral 

nitrogen, the temperature and electron density dependence are given by Page 

+ 1 8 

et al. as; 

7 = K (29) 

The proportionality constant K was obtained from the tabulations of Wilson 

9 

and Hicolet. The computation began by picking first-estimate values of T, 
n^, and n^. These three parameters were evaluated at the same reflected 
shock Hugoniot state (ref. Table A2) and so, for iterative purposes, they 
can be considered as a single variable. For Lorentzian lines and neglecting 
shifts, the spectral absorption coefficient is given by 

* 

The effects of Doppler broadening were negligible for the conditions of this 
experiment . 
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where ^ is the integrated line absorption (ref. Eq. 23 ). 


(30) 


E’ow to generate a spectral intensity^ the continuum must he included. This 
was done by estimating the relative intensity of the continuum^ defined 

as C = Rewriting this expression in terms of absorption coefficients 

(ref. Eq. 19 ) yields 

1 - exp(-iJ,^ 6) = C [1 - exp [- (p,^^ + 6]] (3l) 

Q 

which can be solved explicitly for iJi^. This value ^ together with the line 
absorption coefficient (Eq. 30 ) permits computation of the relative spectral 
intensity profile by the relation 

1 - exp[-(ia.^ + |j,^) 6] 

^Xo 1 - exp[-(iA^^ + 6] 

To compare this theoretical profile to the experimental relative exposure 
profile^ the instrument function cannot be neglected since spectral shape 
is now an important parameter. The above theoretical intensity profile was 
numerically broadened using Eq. (k) and the instrument function shown on Fig. 

1. Thus a theoretical relative exposure profile was generated based on 
assumed values of the parameters T^ n^^ n^^ and C. 

The criterion used to determine the quality of the fit was based on the EMS 
value between about 10-20 points taken from the experimental profile and the 
theoretical values generated at corresponding wavelengths (ref. Eq. 37 on p. 
53 )- The iterative procedure consisted of varying the two input parameters 
(i.e. (t^ n^^ n^) and C) until a minimum RMS value was found. As shown later, 
these minimum values corresponded to an average difference of .01-. 04 units 
of relative exposure which is in good agreement with the uncertainty given 
by Eq. ( 28 ). 
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The dependence of this fitting procedure on such intensity- governing para- 
meters as the nitrogen number density and f-number is rather weak. This can 
be demonstrated readily by again making the assumption that the radiation 
was optically thin. By working backwards through Eqs. (32), (3l), and (30), 
the following relation is readily obtained: 



1 + C(- 


X-Xo 
7 ■ 


2 


1 + (• 


X-Xo 

7 


2 

) 


(33) 


Here it is clear that the optically- thin shape is dependent only on the con- 
tinuum level (which is perturhed)^ the wavelengths (which are considered 
precisely known) ^ and the half half-width which^ from Eq. ( 29 )^ depends 
weakly on temperature and linearly on electron density* Thus the other 
thermodynamic parameters enter only through consideration of a finite opacity 
and this was usually a relatively unimportant (but nonnegligible) effect. 


This shape fitting procedure was applied not only to neutral nitrogen but 
also to the Ea line (6563 A) which was present as an impurity. This pre- 
sented three complications. First the shape is not Lorentzian and is given 
in tabular form by Griem.^^ This was handled by a program option in the 
computer and presented no great difficulty. Secondly^ because hydrogen was 
an impurity^ the number density could not be obtained as before. Instead it 
was determined from a peak spectral intensity obtained using the spectrographic 
results to scale a 6488 A polychromatic measurement of a nitrogen multiplet. 
This was an uncertain number but recall that the number density is relatively 
unimportant provided opacity effects are small. Finally a controversy exists 
as to whether the Griem*^ shapes correctly represent Hct:* Morris et al. 
have calibrated Ecu against the better-known values of Hp and report that 

16 

the Griem values for the Ea half- width are too narrow by 2&fo at n =10 

17 - ^ ^ 

and by l4^ at n^ = 10 cm and cite corroborating results by Wiese at UBS. 

The Eoi shape fitting was done using both shape values but as will be seen, 

this experiment cannot shed light on this controversy. 
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Error Analysis of Line- Shape Thermometry 


For the purpose of this analysis, the measured quantity will he taken to he 
the line half half-width as given hy Eq. (29) • Substituting this expression 
into Eq. (l3) and performing the required differentiation yields: 


T 


d7 


dn 

;i dP 

dK 
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n 
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dp 


K 
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( 34 ) 


Here^ as before^ the temperature sensitivity (irops sharply ^ith temperature 
(ref. Fig. 6) which in this case reflects the fact that the electron density 
is a slowly changing function of temperature at the higher temperatures. The 
proportionality constant K is known to about + 20^ for visible nitrogen 
lines and a + 10^ pressure uncertainty only contributes about + yjo to the 
temperature uncertainty because the derivative is about 0.5 (ref. Fig. 6). 


To obtain the uncertainty in the half half-width measurement^ it is necessary 
to again assume optically- thin conditions. Then Eq. (33) is applicable and, 
by considering 

differential of 7 can be obtained by straightforward differentiation. The 
uncertainty in the relative intensity is given by Eq. ( 28 ) which was shown 
to yield reasonable values and, because the continuum was perturbed during 
the fitting process, it can be considered known to + O.OO 5 units of relative 
intensity. Inserting these uncertainties into the total differential of 
yields ; 


and C as the independent variables, the total 


d7 1 ± O.Oh (l-C) + 0.005 
72 I. 

(^-C) 


dX 

l-lo 


(35) 


This yields the uncertainty in 7 resulting from a single measurement point. 
However the measurement was not done at a single point but rather at many 
throughout the profile. It is assumed that the overall uncertainty is typical 



of the value at the half- intensity point of the line contribution, i.e. 

I. /l^ = C + 0.5 (l“C) where 7 = + (X-X ). With these assumptions, the 

K A.O O 

above equation becomes 


d 7 + 0.04 (l-C) ± 0.005 dX 

7 l-C - 7 ^ 

The wavelength uncertainty will now be neglected; it was carried this far to 
indicate that very narrow lines can cause difficulty. The remaining expression 
yields errors of + at C = 0 ^ + 55 ^ at C = O.5, and + 6 fo at,C = O.75. This 

latter value can be unrealistic because high continuum levels imply large 
overlap and often increase the difficulty of separating one line from another. 
This source of uncertainty was not included in the error analysis. With this 
reservation in mind^ the uncertainty in the half half- width measurement will 
be taken to be + 5^. 


Thus the temperature uncertainty resulting from the shape- fitting process on 
nitrogen lines was + at 12,000 K and this rises to + 7^ at 15,000 K. As 

will be shorn in the next chapter, the broadening and attendant overlap at the 
higher temperatures made it difficult to maintain the + 5^ measurement un- 
certainty. Thus both quantitatively and in the difficulties experienced at 
the higher temperatures, this thermometric technique was similar to that of 
the line intensity measurements discussed previously. 

For Ha, the shapes are somewhat better known but since 7 0'(n^)^^^, the 
temperature sensitivity is reduced by 2/3 with the result that the temperature 
uncertainties from the Ha measurement are quite comparable to those of the 
nitrogen lines discussed above. 

This discussion completes the derivation of the requisite analytical equations 
to measure the plasma temperature and analyze the resultant errors. 



SPECTRAL SURWS 


The next chapter 'will discuss the results of the plasma temperature measure- 
ments. However before the equations of the preceding chapter could be applied^ 
the features of the emergent radiation were identified and it was shown that 
this radiation was indeed constant with time. 

Spectral Identification 

The first stage of the spectral work involved the identification of the 
various spectral features using time-resolved spectrograms. The microdensito- 
meter traces from three such spectrograms are shown on Fig. 7 a<nd are repre- 
sentative of a very low^ an intermediate^ and a high incident shock velocity. 

The ordinate on these traces is optical density and recall that this is 
related to the relative spectral exposure (or intensity) in a nonlinear way 
through the H-D curves. These spectrograms were taken early in the contractual 
year and used glass optics and an unblocked second order. Spectrograms obtained 
later with quartz optics and blocked second order were essentially identical in 
the region shown. The 2500-3^00 A region available with quartz optics showed 
very little of interest - no prominent lines and a smooth continuum. The light 
leak apparent at the red end of the middle trace on Fig. 7 "w-as traced to the 
spectrograph housing and was repaired. 

The three traces on Fig. 7 contain a great deal of qualitative information 
about the plasma radiation. First note that while the temperature increases 
relatively slowly^ the electron density increases by factors of 5 6 between 

each trace. The bottom trace shows a considerable amount of molecular band 
radiation which was primarily from the (1st neg. ) system. The band heads 

of the (2nd pos.) system are also discernable but considerably less pro- 
nounced. The top two traces clearly show the presence of continuum radiation 
but even qualitative statements about the spectral behavior are difficult to 
make - note the similarity between these traces and the spectral sensitivity 
of the 2475 film (Fig. 4). 
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FIGUBE 7 . MICRODENSITOMETER TRACES OF REPRESENTATIVE SPECTROGRAMS 


The character of the line radiation is most interesting. On the bottom trace, 
the lines are generally prominent and narrow. On the middle trace they are 
still prominent but considerably broader and many more lines and multiplets 
are apparent, especially in the 5-6,000 A region. Finally on the top trace 
the lines have become broad, the continuum has come up rapidly, and many 
lines are difficult to see above the continuum. Much of the structure in 
the 5-6,000 A region has disappeared. The 01 (6l57) multiplet is interest- 

ing because it undergoes a 20 A shift toward the red and becomes noticeably 
asymmetric. 

This rapid rise of the continuum as the temperature increases can be qualita- 
tively explained from theory. In the last chapter assuming (l) optically- thin, 
( 2 ) visible region, and ( 3 ) temperatures about 12,000 K, it was shown that the 
temperature sensitivity for the integrated line intensity was 10 (ref. Eq. 27 )* 
But because the lines are becoming broader at the same time, the spectral 
intensity at the line centers does not increase nearly as rapidly. In fact 
for optically- thin Lorentzian lines, the temperature sensitivity of the 
spectral intensity at the line center is only 3*5 which is the difference 
between that of the integrated intensity (lO typ. ) and the half half-width 
( 6.5 typ.). (This can be readily shown from Eq. 30 et al. ) Retaining the 
above assumptions it can be shown that the temperature sensitivity for the 
continuum (both recombination and negative ion attachment processes were 
included) is 10 at 12,000 K. Therefore on a spectral basis as represented 

3.5 

by Fig. 7, the intensity at the line centers increases by T while the 
continuum increases by T and thus the apparent submergence of the lines is 
not unexpected. 

A listing of all the prominent multiplets and many of the lesser ones that 
have been identified at one time or another is given in Table 1. The identi- 
fication process was greatly facilitated by a computer listing of the prominent 
lines and multiplets which gave quantitative criteria not found in finding 
lists or even multiplet tables. This listing was based on the EBS' log gf data 
(if available) together with the half- widths and remaining gf data of Armstrong 
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TABLE 1 . LISTING MD IDENTIFICATION OF PROMINENT AIR MOLTIPLETS 


Radiating 

Specie 

Multiplet 

Number 

Multiplet 

Wavelength 

Theoretical 

Relative 

Strength 

Identified? 

Comments 

NI 

10 

ii-107 

1000 

Yes 

Much broader than predicted 

NI 

5 

4222 

600 

Yes 


01 

* 

6242 

300 

No 

Never sav this one. NBS misprint? 

NI 


4256 

300 

Yes 


NI 

6 

4l46’ 

70 

Yes 


01 

10 

6157 

ko. 

Yes 

Shifts to red, stands out 

NI 

9 

4928 

20 

Yes 

Two lines, 20 A separation 

NI 

21 

6490 

20(x) 

Yes 

Looks broad, many lines 

01 

6 

3692 

20 (x) 

Yes 

Narrow and weak, often buried 

01 

3 

3947 

20 

Yes 


NI 

13 

5349 

20 

Yes 

Many lines, widely spaced 

NI 

Ik 

5295 

20 

Yes 


NI 

20 

6644 

10 

Yes 


01 

5 

4368 

10 

Yes 

Single line 

NI 

* 

5408 

10 

Yes 


NI 

30 

6756 

10(x) 

Yes 

Broad 

NI 

29 

6955 

7 

Yes 

Film sensitivity down 

01 

9 

6455 

7 

Yes 

Broad 

NI 

* 

5177 

6 

Yes 


NI 


4389 

5 

Yes 


NI 

16 

6oo4 

4(x) 

Yes 

Broad 

NI 

22 

6465 

4(x) 

Yes 

21 

Combined with NI above 

01 

65 

6654 

4 

? 

20 

Under NI 

NI 

31 

6723 

3(x) 

9 

Under 

NI 

2k 

5620 

3(x) 

Yes 

Very broad 

NI 

* 

5201 

2 

Yes 


01 

12 

5330 

1 

Yes 

Broad 

NI 

32 

5836 

1 

Yes 


01 

11 

5436 

1 

Yes 



7 

* Tabulated by NBS , not listed by Moore 
(x) Not tabulated by NBS, listed by Moore 
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et al. as tabulated by Wilson and Nicolet. The multiplets were listed 

in order of decreasing multiplet spectral emission coefficient for the 

reflected shock states of Table A2. A spectral ordering was selected because 

of its applicability to film exposure (ref. Eq. 5) although this ordering 

involves the theoretical half- widths which may be uncertain. The ranking 

of multiplets rather than individual lines was a numerical convenience and 

can make a multiplet consisting of narrow^ widely- spaced lines appear too 

high on the list. However, the emphasis here was not on a precise ordering 

but rather for a qualitative list as an aid to identification. While Table 

1 was prepared from the listing at U = 6.8 mm/jisec, the ordering did not 

s 

change much at other conditions and can be considered applicable over the 
range of conditions shown on Fig. 7* Only multiplet wavelengths in the 
2500-7000 A region were tabulated. Most of the IR lines would lie at or 
near the top of the list and hence a future effort will be directed toward 
photographing the IR spectrum. 

It is evident from the comments on Table 1 that all the listed multiplets 
were identified with but three exceptions. Two of these underlie stronger 
multiplets and hence are questionable while the third - the supposedly 
prominent multiplet 01 (6242) - was never seen. While this table is not 

intended to be a complete listing, it includes most of the prominent multiplets. 
There are four other multiplets at wavelengths about 4670, 4485, 4435, and 4350 A 
that, until recently, defied identification- However a comparison with the 
spectral scans by Morris et al. of a nitrogen arc definitely establishes 
these to be nitrogen. 

There are several mis- identifications on Fig. 7 that have been rectified in 

10 ^ 

Table 1. For example, on Fig. 7^ "the 01 (39^7) ought to be 01'^, the 4670 A 

17 l8 

identification of 01 and 01 is erroneous - the lines are far too narrow 

2 30 
(also see above), and the 01 identification at 6750 ought to be . 

The only impurity lines ever positively identified in emission were those of 
hydrogen - Hh ( 6563 ) was ever-present and prominent and Hp (486l) was 
apparent on most spectrograms. The hydrogen probably came from water vapor 
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either present on the shock tube walls and not removed by pumping or present 
in the test gas where it was ineffectively trapped. The shock tube was 
typically evacuated to below 0.001 torr while the initial pressure was 0.2 
torr. The niJunber densities obtained during the Ha fitting procedure dis- 
cussed earlier yielded typical hydrogen concentrations of a small fraction 
of a percent. 

If other impurities were present^ their radiative signatures were always 
buried in the noise level of the film. This was an unexpected result since 
no claim is made for extreme shock tube cleanliness although the tube was 
cleaned by dry- wiping between each shot. This result means that the spectrally- 
integrated radiation was not significantly affected by impurities, at least in 
the 25 OO- 7 OOO A spectral interval. 

On one shot the rapid-closing shutter failed to close and the spectrogram 
was exposed by the entire radiative history of the reflected shock region- 
The resultant overexposed film revealed a number of absorption lines and the 
lines of Fe, Cr, Ni, and A1 were positively identified- The first three 
probably came from the stainless steel shock tube walls while the only 
aluminum known to be in the system was that of the trigger wire in the 
driver. 

Another purpose of the spectral surveys was to locate intervals which were 
free of line radiation and hence could be used to measure the continuum 
intensity via the polychromator- Unfortunately, as a scrutiny of Fig. 7 
will show, there was no such interval where one could unequivocally state 
that an intensity measurement was that of only the continuum- Nevertheless, 
since two extra polychromator channels were available, exit slits were 
placed near the mercury green line (5^60 A) where one might hope to measure 
the continuum and at 7032 A which was a fool’s paradise since the film sensi- 
tivity (Fig. 4) was too far down to see anything. These two channels are 
called "continuum?" because about the best that can be said is that they 
yield an upper bound. Their main function was to provide a basis for 
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comparison against the continuum intensity levels that were a hy-product of 
the line intensity measurements. 

Proof of Steady State 

This section is devoted to a demonstration of the temporal behavior of the 
traces from the photomultiplier tubes in the polychromator and spectrograph 
(zero-order). The object is to show that the radiation was fairly constant 
with time and moreover that the behavior of the various spectral features 
was quite similar. 

This latter fact is shown by Fig. 8 which contains traces from an early 
polychromator configuration where the exit slits were located to obtain the 
time behavior of six different spectral features - molecular radiation, lines 
of and 01, the continuum?, an impurity line, and the zero-order spectrum. 
All six traces are remarkably similar in their temporal behavior and this 
was taken to indicate the presence of chemical equilibrium. The fact that 
the nonspectral zero-order trace was among these six means that a white- 
light photomultiplier trace is a valid measure of temporal variations in 
plasmas of this type. The 30-^0^ droop began at about 30 \xsec which was 
the end of spectrograph exposure time. The trace shows a small signal 
about 10 |j.sec prior to the arrival of the reflected shock which corresponds 
to the passage of the incident shock through the field of view. The traces 
shown on Fig. 8 are typical of a series of four shots at various shock 
velocities taken early in the experiment to establish the presence of 
chemical equilibrium and show that time- integrated spectroscopic recording 
of film was a valid technique. 

To further demonstrate the temporal behavior, various photoelectric traces 
from twelve recent shots are shown on Fig. 9 where they are arranged in order 
of increasing incident shock velocity. The traces are principally from either 
the 5^69 A continuum channel of the polychromator or the model profile which 
was the white- light monitor taken between the energy detectors (see ref. 1 for 
more details). If neither of these were available, then the zero-order trace 
of the spectrograph was used - this is less desirable here because the rapid- 
closing shutter temiinates the trace at the end of the exposure time. When 
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27 A bandpass 



Fig- 8 Photomultiplier Traces from Polychromator on Shot #590 
(U« = 6.85 mm/|jLsec) 
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SHOT 603 
(5469 A) 
6.95 


SHOT 601 


SHOT 602 
(5469 A) 
7.50 ^ 




SHOT 578 

SHOT 605 / 


(O-order) 

(5469 A) / 

\ 

\ 

7.74 


SHOT 585 
(mod. pro.) 
7.42 i 


SHOT 606 
(5469 A) 


SHOT 594 
(mod. pro.) 
7.42 


SHOT 604 
(O -order) 
7.44 


SHOT 586 
(mod. pro.) 
8.30 . 


SHOT 597 
(mod. pro.^ 
8.8 


5 (jisec (all sweeps) 


Dashed curves are zero-order traces from spectrograph 
Incident shock velocities are in mm/(xsec 


Fig. 9 Various Photoelectric Traces from Twelve Shots 






available, the zero- order traces are superimposed on the others to indicate 
these times. Because of the similarity between the various photoelectric 
traces on any one shot, the choice of the particular one used here is un- 
important . 

Although shot-to-shot variations are clearly apparent on Fig. the broad 

trend of these traces shows that they are generally flat for U < 7-5 ™i/ 

s 

p,sec while those at higher incident shock velocities generally show a gradual 
decrease with time. This may very well have been caused by radiative cooling 
although, because the radiation is such a sensitive function of temperature, 
the resultant temperature decrease is quite small. For example, the tempera- 
ture sensitivity of the visible continuum at l 4,000 K (U = 8.3 nim/|Jisec) is 
7.7* Assuming a liberal decrease in the 5^69 A continuum trace over 
25 iJLsec, the corresponding drop in temperature would be only 15/7»7 or 2^. 
Further because the thermometric measurement techniques discussed earlier 
were both time- integrated to a certain extent, the measured temperatures 
represent a time- averaged value which would deviate from the original value 
by about 1 ^. Considering that the precision of the measurements themselves 
was about ± 4 ^, this variation is insignificant. Thus for the purpose of 
measuring plasma temperatures, the assumptions of chemical equilibrium and 
steady state conditions were valid. 

The test times behind the incident shock were determined through other photo- 
electric monitors (see ref. l) and were generally 20-30 jjbsec. For these 
experiments, the geometry was such that these times are very nearly equal to 
the times required for a sonic disturbance from the reflected shock- contact 
surface interaction to propagate into the field of view. These would be 
minimum values for the reflected shock test times and are compatible with the 
25-30 M^sec spectro graphic exposure times. 

Selection of Plasma Thermometers 

On the basis of scans like those of Fig. 7, the various lines and multiplets 
used to measure the plasma temperature both through shape and intensity 
measurements were selected. 
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For intensity measurements^ well-known^ prominent atomic lines or multiplets 
were desired. The (i )-935 A) line and the 01 ^® (6157) and ( 6490 ) 

multiplets were selected since they possessed the fewest disadvantages. The 
(^935) line was bright and narrow^ but was located in a spectral region 
where film sensitivity variations could be important (ref. Fig. 4 ). The 

prominence of the 01 ^^ ( 6 l 57 ) overshadowed its asymmetry and red shift. The 

21 T 

WI ( 6490 ) multiplet was not always prominent and is not listed by KBS^ 

but it was close to Hd where an absolute intensity was needed and it was 

tabulated by Wilson and Mcolet.^ In retrospect the ( 6644 ) multiplet 

21 

might have been a better choice since it was listed by but WI was 

somewhat more easily separated from other lines. The polychromator exit 
slit widths were selected to be about six times the line or multiplet half 
half-width at nominal conditions and are tabulated below. The wavelengths 
are those of the centerlines and AX is the full width of each exit slit. 


Spectral feature 

\,A 

AX, A 

9 

WI^ 

4935 

4.6 

Continuum? 

5469 

16.6 


6163 

45.8 

H 

C\J 

H 

6488 

47.1 

Continuum? 

7032 

16.7 


For shapes^ the candidate had to be a single line in addition to being well- 
known and prominent. WI^ ( 4935 ) ^i^et this added restriction and Hd (6563) 
was considered well-known despite the shape controversy mentioned earlier. 

The better- known HP ( 486 l) was not prominent and could only be used on 
one shot when an unusually high (l.3^) hydrogen concentration made it usable. 
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RESULTS 


The results of the plasma temperature measurements are presented and dis- 
cussed in this chapter. They were obtained on the six shots for. which both 
spectrographic and polychromatic results were available. These shots had 
incident shock velocities in the range of mm/p»sec and yielded average 
measured temperatures from 10,750 to 13^000 K. Examples of the techniques 
and interesting qualitative results will be presented first. 

Examples of Shape-Fitting 

In this section are presented examples of the shape- fitting process taken 
from these six shots. The lines considered were the 4935 A line of 
Ha (6563 a) and, on one occasion, Hp (486l A). The criterion used to 
determine the best fit between the computer- generated shape and the experi- 
mental profile was that the fit criterion 6e be a minimum where; 

6e = ^r 2 [e(\ ) - e(X ) f\ (37) 

^ i=l expt. ^ theor. 

Here e(X^) refers to the experimental and theoretical relative exposure 
values at a wavelength The minimum value of the fit criterion is shown 

on the examples which follow. 

Figure 10 shows the measured profile and the best- fit theoretical data points 

for Ha on one of the slower incident shock velocities using the shape taken 

from Griem. The minimum value 6e = 0.033 "^^as higher than most for this 

experiment. The fit using the Morris et al.^^ correction to the Griem half 

half- width (the width at every point was changed by this amount) would be 

similar but the resultant electron density would be lower. The results for 
9 

NI on this same shot are given in Fig. 11 which also contains the instrument 
function to indicate that it must be included in the data reduction procedure 
for shape measurements on narrow lines. As indicated by the value 6e = 0.011 
the theoretical points fit the measured profile quite well. 
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RELATIVE SPECTRAL EXPOSURE 



0 +4 +8 +12 


A - Aq , A (A^ = 6563 A) 

Fig. 10 Example of Shape-Fitting Process for Ha on Shot #601 
(U„ =7-00 mm/(isec) 






Following this shot (which was the first of the six) it was decided to fit 
the entire profile rather than just half. To fully utilize the features in- 
herent in the RMS fit criterion (Eq. 37)^ the measured profiles were read 
directly into the computer as taken from the microdensitometer traces even 

though these profiles contained obvious noise. An example of a noisy profile 

9 

is shown by Fig. 12 for at an intermediate shock velocity. Here the 

value 6e = 0.039 rather high which was to be expected from the nature 

of the measured profile, but the resultant electron density yielded a 

temperature in good agreement with the other measurements as will be shown 

later. The exit slit width of the polychromator AX is shown in this figure 

for informative purposes only - the polychromatic result was not used in the 

9 

fitting process for NI and only in a relatively unimportant way for Ho'. 

9 

The NI profile for the highest incident shock velocity of these six shots 

is shown by Fig. 13 which also indicates the rapid rise of the underlying 

continuum. The attendant high electron density broadened the line to the 

point where the wing of the ^915 A line in the NI^ multiplet overlapped 

noticeably. This overlapping was accounted for (approximately) by lowering 

the profile on the blue wing as shown. Another and probably better approach 

9 

would be to fit both lines of the NI multiplet simultaneously. This 

technique would be applicable here because there are no other adjacent lines 

to cause difficulty (NI^ is singular in this respect). Even though the 

continuum level was relatively high, the fit was quite good as indicated 

by the low value of the fit criterion (6e = 0.006). These last three figures 

9 

also indicate that NI does indeed have a Lorentzian line shape. 

As mentioned earlier, HP was often only barely discernable above the con- 
tinuum. This happens because it is inherently a weaker line than Hu (the 
maximum absorption coefficient is I/30 of that for Hu) and is located in a 
region where the air continuum is higher. However, on one shot the H^ con- 
centration was abnormally high 1-3^) and HP stood well above the con- 
tinuum. This profile and the results of the shape- fitting process are shown 
on Fig. lU. The calculations were done by hand using the shape given by 
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Grlem^*^ and the effects of the instrument function and opacity were neglected. 
The ahcissa of Fig. l4 indicates that the instrument function was relatively 
quite narrow while an intensity calculation revealed that the peak emissivity 
was about 0.10 so opacity effects were small. 

This figure shows that the best fit profile corresponded to an electron 
density of 5 x 10"^ Portions of the shapes for electron densities 

perturbed by + 20^ are shown to indicate the sensitivity of the fit. 
is a well- studied line whose shape is reported by Griem to be good to 3^ 
for n^ = 3 X 10^^ - 3 x 10^"^ cm ^ and T = 10-20^000 K. Because of this 
low uncertainty on the theoretical shape^ an error analysis similar to that 
done earlier indicates that HP was a + 2^ thermometer on this shot. This 
suggests another possible thermometric technique that will be further discussed 
later. 

The reason why the Ecx and HP shapes cannot be used simultaneously is 

21 

typified by this shot. The relative exposure profiles for NI and Ha 

revealed that the Ha peak intensity was 3*7 times that of the fairly-bright 
21 

HI multiplet. Scaling up the polychromatic result gave an emissivity very 
near unity. With this high opacity^ the shape-fitting process becomes far 
too dependent on an intensity measurement^ e.g. lowering the Ha peak 
intensity by 20^ doubled the resultant best- fit electron density. This is 
obviously an unsatisfactory situation and hence the Ha shape was essentially 
useless on this shot. This effect could be mitigated by shortening the path- 
length and later this approach will be further discussed. However^ for the 
30.5 cm pathlengths of this experiment^ it was impossible to measure both 
Ha and HP simultaneously under conditions where the shape- fitting process 
was meaningful for both lines. 

Examples of Line Overlap 

9 

Figure I 3 revealed that a small overlap occurred between the lines of HI at 

the higher electron densities. While this did not cause much difficulty on 

9 

HI ^ line overlap seriously limited the usefulness of the polychromatic results 
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for and 01^^. This is illustrated hy Fig. 15 which contains the rela- 

tive exposure profiles for the 6^4-10-6570 A region for three shots - an inter- 
mediate and the two extremes of incident shock velocity. This spectral 
region contains one 01 and three ITI multiplets and Hor. The wavelengths for 
these are indicated at the bottom of the figure as is the position and width 
of the appropriate polychromator exit slit. By superposing the three profiles^ 
the effects of line broadening^ shifting^ and overlapping and the rising of 
the continuum become most apparent. 

. 9 21 

It is evident that only 01 ^ NI , and Ea are the important radiators in 
this region^ but despite this simplification it becomes increasingly difficult 
to separate the contributions from these three multiplets at the higher in- 
cident shock velocities. In fact, no attempt was made to separate them on the 

upper profile (shot 6o6) because the errors on the resultant temperature ob- 
21 

tained from the NI integrated intensity would be unreasonably large. This 

typifies the difficulties encountered in attempting to measure the intensity 

of a line or multiplet that is not well isolated. Figure 15 also typifies 

the difficulties in attempting to fit the shape of Ho' since the profile is 

21 20 

strongly influenced by the overlap from NI (and from NI on the red side 
of Hh). 

Examples of Combined Results 

Examples of the combined results from the polychromator and the spectrograph 
are given by the following three figures (Figs. l6-l8) which again represent 
an intemediate and the extremes of incident shock velocity. The solid 
symbols on these figures represent the polychromatic results (ref. Eq. lO) 
for the three line and two continuum measurements using the exit slits 
tabulated on p. 52. As discussed earlier, these solid symbols represent a 
spectral intensity averaged over the spectral width of the exit slit. 

For each of the three polychromatic measurements which included a line or 
multiplet, combination with the relative exposure profile from the spectro- 
graph yielded a peak spectral intensity (Eq. 12) if the effect of the instrument 
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Fig. l6 Combined Results for Shot ifGOS (Uo = 6.95 mm/iisec) 
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Fig. 17 Combined Results for Shot #605 (U = 7.36 mm/nsec) 
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Fig. l8 Combined Results for Shot #6o6 (U_ = 7-99 mm/lisec) 




function was neglected. Also the shape- fitting process yielded a value for 
the continuum level (ref. p. 38 )- Both of these values are indicated by the 
open symbols on Figs. 16-I8 and represent true spectral intensities. 

The main purpose of these plots was to provide a graphical comparison between 
the inferred continuum level from the three line intensity measurements and 
the two direct measurements of the continuum. This provided a consistency 
check and on one occasion uncovered an error in the microdensitometric process. 
For reasons discussed in the next section^ the continuum prediction was fit on 
the 4935 A continuum intensity and the resultant curves are shown. These 
curves were generated from a smeared continuum model where the photoionization 
edges are lumped and hence detailed spectral features disappear. Hence little 
more can be said about this curve except that the other continuum points lie 
quite close or slightly above it. The important fact is that these points 
showed consistency from shot to shot. 

The Planck function^ evaluated at the temperatures associated with the con- 
tinuum predictions^ is also shown on these three figures. It is noted that 
the Planck function changes relatively slowly with temperature. The much 
more rapid rise of the continuum and the attendant approach to the Planck 

function is shown quantitatively by these three figures. On Fig. l8^ the 

9 21 
NI peak emissivity was 0.29 while the NI value was about 0.5. Obviously^ 

these peak emissivities can be reduced by shortening the pathlength. 

However of much more importance to the validity of the temperature measure- 
ments is the fact that the continuum intensity is approaching or even exceed- 
ing the peak intensities of the lines. This is demonstrated by mentally 
picturing where a curve of twice the continuum prediction would fall on 
Fig. 18. How since the polychromatic measurement includes contributions 
from both the lines and the continuum^ high continuum levels mean that the 
line intensity is becoming the difference between two large numbers. As 
pointed out in the error analysis, this leads to large uncertainties in the 
integrated line intensity and hence large errors in the resultant plasma 
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temperature. It is the rapidly- rising continuum^ and not the increased 
opacity, that limited the usefulness of this thermometric technique. 

The Thermometric Use of the Continuum 

A possible thermometric technique that has not been discussed as yet involves 
the continuum itself. This is basically attractive because of a high tempera- 
ture sensitivity of the continuum and it is this that has made the line 
measurements increasingly unattractive. However, the theoretical continuum 
cross sections are generally quite uncertain and limit the usefulness of 
such a thermometer. 

l4 

But Morris et al. have extensively studied the continuum of nitrogen at 
4955 A and recommend that their values of absolute intensity be used to deter- 
mine a plasma temperature. Accordingly, the values for the spectral continuum 
intensity of air at 4955 A taken from the theoretical predictions (see ref. l) 
which used a smeared continuum model were compared against the detailed results 

l4 

of Morris et al. for nitrogen at corresponding temperatures and pressures. 

The agreement was remarkable (and very probably fortuitous). The predictions 
were about 10^ high in the 10,500-12,000 K range (which is to be expected 
when comparing air to nitrogen) and essentially in agreement in the 13,000- 
l4,000 K range of reflected shock conditions- 

The result of the integrated intensity measurement and shape- fitting process 

on (4935 a) gave, as a byproduct, a value for the continuum spectral 

intensity. Also the continuum predictions, when fit on these values, were 

consistent with the other four continuum measurements as shown by Figs. 16- 

18 and thus increased our confidence in these 4935 A values. Thus because 

good experimental values for the 4935 A continuum spectral intensities were 

available and because the predictions were essentially in agreement with the 

l4 

careful experimental studies of Morris et al., it was decided to evaluate 
a plasma temperature based on these predictions. The results will be dis- 
cussed in the next section. 
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Although placing an exit slit over a line is admittedly not an optimum way 
to obtain the continuum intensity^ this thermometric technique was applied 
a posteriori and improvements in technique for future measurements could 
obviously be made. The large temperature sensitivity inherent in this 
method is demonstrated by the large difference between the two continuim 
predictions placed on Fig. 17 . 

Results of the Plasma Temperature Measurements 

The culmination of the spectroscopic effort is represented by Fig. 19 which 
shows the results of the various plasma temperature measurements plotted 
against the gasdynamic value obtained from the incident shock velocity and 
initial pressure (ref. Table A2). To aid in the assessment of the data 
scatter^ + error bars were drawn about the line representing temperature 
equality. These bars represent a reasonable estimate of the experimental 
errors for the temperature range shown and, with the exception of shot 602 , 
the scatter between all the results on any one shot lie within the heights 
represented by these bars. 

The important conclusions to be drawn from this figure are first that the 
measured temperatures do not cluster about the gasdynamic value as would be 
expected if the gasdynamic temperatures were representative of the state of 
gas and if the measurement errors were randomly distributed. Also some 
erratic behavior is indicated - e.g. the mean of the measurements on shots 
601 and 605 are nearly the same while the gasdynamic values differ by almost 
1,000 K. This erratic behavior also suggests that the gasdynamic value does 
not adequately define the state of the test gas. 

Figure 19 contains all the results that could be derived from the six shots 

for which both polychromatic and spectrograph! c results were available. 

Since for each shot there were 3 line intensity measurements (RI^, 01^^, 

21 9 

WI ) and 2 shape measurements (NI^ and Ha), this should yield 30 data 
points. In actuality there were only 22. The lost points were caused by a 
misaligned exit slit (l), a bad spot in the negative caused by faulty 
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development (l)^ off-scale polychromator traces ( 3 ) caused by an unexpectedly 
hot shot, and the inability to utilize the Hh, and 01^^ data on shot 

606 because of the high continuiam level (ref. Fig, 15)* 

On shot 603 , the unusually- high hydrogen concentration made Ea unusable but 
this in turn yielded a good profile for HP (ref. Fig. l4) which, as dis- 
cussed earlier, was a 2^ thermometer. The resultant temperature is gratify- 

9 

ingly close to the mean of the other measurements. Also the El intensity 
measurements made possible the continuum thermometer discussed earlier. 

These points are included on Fig. 19 and also lie quite near the mean of the 
measurements. Note also the good correspondence between the HP and the 
4935 A continuum results on shot 603 . 

The above statements are clarified by Fig. 20 which compares the deviation 

of the various measurements from the arithmetic mean of all the results 

(including the 4935 A continuum but excluding the gasdynamic value) available 

on each shot. The HP value and the four 4935 A continuum measurements lie 

quite close to the mean value. It is interesting to note that the scatter of 

9 

the HI integrated line intensity measurements is about twice that of the 
4935 A continuum. This must be explained at least in part by a higher tempera- 
ture uncertainty for the continuum measurements since both absolute intensi- 
ties were obtained from the same polychromatic result. 

Another reason for preparing Fig, 20 was to locate possible one-sided devia- 
tions which would point to a systematic error. None is clearly apparent 

although the 01^^ and NI^^ measurements do appear slightly below the mean 

9 

while the NI shapes all lie above the mean. These however are within the 
expected errors from uncertainties in the respective f-numbers and half half- 
widths. 

More important is the last column of Fig. 20 which shows the deviations of 
the gasdynamic values of the plasma temperature. These are definitely one- 
sided and average about 4^ above the mean of the measured temperatures. This 
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DEVIATION FROM MEAN OF MEASURED TEMPERATURES (%) 



603 6.95 10,760 
A 601 7.00 11,160 
O 605 7.36 11,380 
□ 604 7.44 11,385 
O 602 7.50 12,025 
/i 606 7.99 12,950 

Fig. 20 Comparison of the Results from the Various Thermometric Techniques 


T1 


is a significant deviation and indicates that the actual reflected shock 
temperature was somewhat helow the gas dynamic value. The erratic behavior 
of the gasdynamic value is exemplified by the fact that there is no corres- 
pondence between the magnitude of the deviation and the incident shock 
velocity. 

All the Ea results on Fig, 20 and most of those on Fig. 19 were obtained 
using the shapes recommended by Morris et al.^^ As shown on Fig. 19^ the use 
of the Griem^^ shape would raise the measured temperatures about 4^. As 
indicated by Fig. 20^ the results using the Morris shapes were not sufficiently 
one-sided to make one shape a clear choice over the other. In fact^ to center 
the Hcv temperatures about the zero deviation line would require an inter- 
mediate shape about halfway between the two. 


The intermediate shapes recently computed by Kepple and Griem meet this re- 
quirement and probably should be used in the future. The reference is Kepple^ 
Paul C., "improved Stark Profile Calculations for the First Four Members of 
the Hydrogen Lyman and B aimer Series^" University of Maryland Tech. Report 
#831, 1968 . 
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CONCLUSIONS AND RECOMMENDATIONS 


The measured temperatures of this experiment are in reasonable agreement 
■with each other vithin the limitations imposed by the expected errors. Ho-w- 
ever the mean of these measured values lies about belo-w the gasdynamic 
temperature. Further^ although the number of shots is not as high nor the 
scatter on any one shot as low as desired, there are indications that the 
actual temperatures may be randomly scattered in a manner not predictable 
from the incident shock velocity. 

To place the results of this experiment in the proper perspective, recall 
that the temperatures were measured to shed light on the deviations between 
experiment and theory for total radiation measurements. At 12,000 K, which 
is representative of the measurements described herein, the theoretical total 
radiant intensity predictions lie above the experimental results by a factor 
of about 1.5 (ref. Fig. Al). Since here the total intensity varies as the 
tenth power of temperature, this factor of I .5 could be explained by a 
temperature difference. 

Thus the average difference observed between the measured and the gas- 
dynamic temperatures is certainly in the right direction and of about the 
correct magnitude to explain the deviation between theory and the experi- 
mental total radiation measurements. However it is difficult if not impossible 
to justify any stronger conclusions on the basis of six ± 4^ measurements when 
a 4^ effect is all-important. 

The temperature measurements reported herein were not significantly affected 
by radiative cooling phenomena. The analysis on p. 53^ which was based on 
actual photoelectric traces obtained at temperatures higher than those 
measured herein, showed that cooling effects were less than Vjo. 

In this experiment, temperatures were only measured to 13^000 K which fell 
short of the desired goal of l6,000 K. The combination of the following two 
effects limited this range : 
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1. The rapid rise of the continuum made it difficult to perform 
adequately precise measurements of the radiative quantities, 

2. The temperature sensitivity of the various thermometric techniques 
decreased sharply. 

Thus in short it rapidly became more difficult to make even less-precise 
measurements • 

ITevertheless it is recommended that additional temperature measurements be 
made and further that these be as precise as possible. Clearly this will 
require a preliminary analytical effort. However experience in these matters 
has been gained and the chapter on plasma thermometry outlines a systematic 
method of approaching the problem. The following thoughts ought to be kept 
in mind; 

1. At the higher temperatures^ radiative cooling effects may become 
important and hence not all the techniques should require that the 
plasma be in a steady state. (A direct, photoelectric measurement 
of the U 955 A continuum is a possibility here.) 

2. The plasma could be seeded with another specie which would offer 
attractive advantages (e.g. the shape of HP). 

Also the existing experimental techniques should be improved to include a 
measurement of the plasma pressure and the use of a step wedge with finer 
steps to yield more points on the H-D curve at higher relative exposure 
levels. In addition the pathlength should be shortened to minimize opacity 
effects and all spectroscopic apparatus should then be located on the same 
side of the shock tube. 

This report is concluded with a brief comment about the ramifications that 
these improved and extended property measurements would have on the measure- 
ments of the total radiant intensity. The erratic behavior of the gasdynamic 
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temperature encountered during this research strongly suggests that plasma 
property and total intensity measurements -will have to he performed simul- 
taneously. While this would require an additional effort^ the present cavity 
gage configuration (see ref. l) was designed for simultaneous measurements 
and the THT gages are proven and reliable energy detectors. Thus this task 
would not require an inordinate amount of time and would do much toward 
validating the current state-of-the-art predictions for the radiant intensity 
of high-temperature air. 
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APPEKDIX 


This appendix contains two figures and two tables included to supplement the 
discussion in the main text. 

The figures present the salient results of the latest spectrally- integrated 
radiant intensity measurements taken using energy detectors on both this and 
prior NASA contractual efforts. This work is described in ref. 1. Figure 
A1 presents the experimental results and a theoretical prediction for the 
total (i.e. windowless) radiant intensity of air while Fig. A2 presents the 
results and a prediction for the radiation transmitted by a quartz window. 

Table A1 contains numerical values for the total radiant intensity prediction 
at the conditions of this experiment. These predictions are felt to be repre- 
sentative of the current state-of-the-art. The origin of these predictions 
is given in ref. 1. It is noted that the values in the first column (line^ 

< 2;,000 a) do not always vary smoothly with pathlength. These irregularities 
were caused by switching from an analytic expression to an exact numerical 
frequency integration for the equivalent width of certain line groups after 
the lines had become strongly overlapped. This caused 10^ errors in the first 
column and^ of course^ much lesser ones in the total intensity column. 

Finally Table A2 contains a listing of the thermodynamic states and the species 
composition in the reflected shock region as computed from the Pankine-Hugoniot 
equations and an IMSC thermochemical equilibrium program. The data from this 
table were used to prepare the theoretical predictions shown on Figs. A1 and 
A2 and certain parts were used in the plasma thermometry discussed in the main 
text of this report. 
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Fig. A1 Experimental Results and Kieoretical Predictions for the Total 
Radiant Intensity. 
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TABIE A1 

RADIANT INTENSITY PREDICTIONS AT CONDITIONS IN THE REFIECTED SHOCK REGION 


Conditions 

Path 
6, cm 

ALL INTENSITIES HAVE UNITS OF w/ci^-SR 

Line 

<2000A 

Line 

>2000A 

Con- 

tinuxHn 

<2000A 

Con- 

tinuum 

>2000A 

Neg. 

Ion 

Molecular 

Band 

Total 

( 6, 6 imn/ijisec) 

1 

10.1 

7.79 

8.21 

2.29 

3.71 

1.74 

33.9 

9,762 K 

3 

15.8 

17.0 

13.6 

6.81 

11.1 

5.22 

69.6 

2,l84 atm 

5 

19.4 

24.4 

16.5 

11.3 

18.5 

8.70 

98.8 


8 

23.4 

33.4 

19.1 

18.0 

29.7 

13.9 

137 


12 

27.8 

43.6 

21.0 

26.8 

44.5 

20.9 

185 

(6.8 inm/|j,sec) 

1 

33.3 

20.0 

25.2 

8.09 

8.75 

2.04 

97.4 

10,490 K 

3 

51.3 

44.0 

4o.l 

24.0 

26.2 

6.12 

192 

2.380 atm 

5 

62.1 

62.9 

47.6 

39.6 

43.8 

10.2 

266 


8 

73.7 

85.6 

53.9 

62.8 

70.0 

16.3 

362 


12 

88.7 

112 

58.4 

93.1 

105 

24.5 

481 

(7.2 iran/jjLsec) 

1 

172 

73.9 

120 

46.1 

29.0 

2.26 

443 

11,710 K 

3 

248 

163 

179 

135 

86.8 

6.77 

819 

2,791 atm 

5 

304 

229 

207 

220 

145 

11.3 

1,120 


8 

299 

309 

230 

344 

232 

18.1 

1,430 


12 

348 

397 

248 

504 

348 

27.1 

1,870 

(7.8 mm/p,sec) 

1 

746 

237 

484 

215 

85.4 

2.00 

1,770 

13,070 K 

3 

799 

507 

694 

611 

256 

6.00 

2,870 

3.437 atm 

5 

961 

691 

789 

983 

427 

10.0 

3,860 


8 

1, 120 

902 

874 

1,510 

682 

16.0 

5,110 


12 

1,290 

1,110 

961 

2,170 

1,020 

24.0 

6,590 

(8. 4 mm/pisec) 

1 

1,300 

849 

1,210 

568 

167 

1.30 

4,090 

I4,l60 K 

3 

1,890 

995 

1,700 

1,570 

502 

4. 50 

6,650 

4.083 atm 

5 

2,190 

1,300 

1,930 

2,490 

835 

7.50 

8,740 


8 

2,440 

1,600 

2,150 

3,750 

1,330 

12.0 

11,300 


12 

2,780 

1,800 

2,410 

5,280 

2,000 

18.0 

14,300 

(9-0 mm/ij,sec) 

1 

3,110 

770 

2,320 

1,070 

246 


7,520 

15,080 K 

3 

3,390 

1,470 

3,24 o 

2,910 

738 

- 

11, 800 

4.601 atm 

5 

3,910 

1,740 

3,690 

4,530 

1,230 


15,100 


8 

4,220 

2,070 

4,160 

6,690 

1,970 


19, 100 


12 

4,820 

2,330 

4,720 

9,220 

2,950 

- 

24,000 

(9-6 mm/p,sec) 

1 

4,820 

1,010 

3,710 

1,590 

291 


11, 400 

15,870 K 

3 

5,150 

l,84o 

5,220 

4,240 

872 


17,300 

4,889 atm 

5 

5,550 

2,110 

5,970 

6,530 

1,450 


21, 600 


8 

6,220 

2,420 

6,780 

9,510 

2,330 


27, 300 


12 

6,910 

2,630 

7,710 

12,900 

3,490 

- 

33,700 

(10.2 mm/tisec) 

1 

7,230 

1,280 

5,830 

2,300 

341 


17,000 

16,700 K 

3 

7,500 

2,090 

8,250 

6,030 

1,030 

... 

24,900 

5.374 

5 

8,060 

2,460 

9,*^90 

9,160 

1,710 

- 

30,900 


8 

8,94 o 

2 , 7 io 

10,800 

13,200 

2,730 

.. 

38,300 


12 

9,810 

2,870 

12,300 

17,700 

4,090 

- 

46,800 
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TABLE A2 

CALCULATED THEEMODYMMIC STATES OF THE REFLECTED SHOCK REGION 
(P^ = 0.2 Torr; T^ = 29 h Gas ’ J 9 io N^;, 2lfo 0^ by pressure) 


Us 

nmi/ 

^sec 

"5 

“K 

"5 

atm 

^5 

gm/ cc 

U5 

cal/ gm 

NUMBER DENSITIES IN PARTICLES/(CM)^ 

e 

N 

0 

n"^ 

o"^ 

^2 

°2 

n“ 

0“ 

6.6 

9,762 

2.184 

3.917'^ 

1.107"^^ 

2.07'*'^^ 


3.42^^"^ 

1.7>.*18 

3.24"^^^ 

1.43*18 

7.24"^^^ 

1.3 o ""2 

5.41'^^^ 

6.8 

10,490 

2.381 

3.912'^ 

- 

1.181*^ 

3.86"^^^ 

1.25^^® 

3.38’^^'^ 

3.28"^^^ 

5.78'^^^ 

6.3 i ""5 

5.12'^^^ 

2.34"^^^ 

7.90“^^^ 

7.2 

11,710 

2.791 ! 

3.976"^ 

1.331'"^: 

9.24“^^^ 

1.24'*'^® 

3.37"^^ 

7.93'"^^ 

1.33*18 

1.99“^^^ 

3.25^^^ 

5.1 i "12 

QIQ 

7.8 

13,070 

; 

3.437 

4. 144' 5 

1.566"^^ 

2.02’^^'^ 

1.20"^^® 

3.36^"^ 

1.73^^^ 

2.89^^^ 

6.71*1'* 

2.17'^^^ 

9.84’^^ 

2.06"^^3 

8.41 

l4,l6o 

4.083 

4.275"^ 

1.816"^^ 

1 1 

3.34"^^' 

1.13*“ 

3.28*11 

2.86^^^ ! 

4.86''^^ 

3.00'^^^ ; 

1.57^^^ 

1.4 i "13 

2.64^^3 

9.0 

15,080 

4.601 : 

4. 248' 5 

> 2.078"^^: 

4.69""^ 

1.01*18 

I 

3.04'^^'^ 

3.99''^'^ 

7.03^^^ i 

1.44^^^: 

1.10"12 

1.65"'-3 

2.87'^^3 

9.6: 

15,870 

4.889 

4.028"^ 

2.350^^ 

■ 5.83^^^ i 

8.43'"^^ 

2.64'^^’^ 

4.92’'^^ 

9.07"^^^ ' 

6.78'^^^ i 

7.03"^^^ , 

1.6 i "13 ^ 

2.69'*'^^ 

10. 2 1 

16,700 

1 

5.37^ 

3.9^7‘^ 

2.646^^! 

. 7.18’^^'^ 

1 1 

7.08'^^'^ 

2.32'^^'^ 

6.02'^^'^ 

1.16*11 ; 

3.27"^^^ 

^.57"^/ 

1.56"^3’ 

2.52‘^^3 

10.8 

;17,610 

6.0371 

3.946"^ 

2.967"^^ 

8.74^^^ 

i 

5.83"^^ 

1.99'^^'^ 

7.26'*'^'^ 

i 1.48^1^ 1 

1.53'"^^ 
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